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We consider a Markovian jumping process with two absorbing barriers, 
for which the waiting-time distribution involves a position-dependent co- 
efficient. We solve the Fokker-Planck equation with boundary conditions 
and calculate the mean first passage time (MFPT) which appears always 
finite, also for the subdiffusive case. Then, for the case of the jumping-size 
distribution in form of the Levy distribution, we determine the probabil- 
ity density distributions and MFPT by means of numerical simulations. 
Dependence of the results on process parameters, as well as on the Levy 
distribution width, is discussed. 

1. Introduction 

Transport processes in physical systems are usually considered in the 
diffusion limit of large distance at long time, fn realistic situations, how- 
ever, the available space is finite and this must be taken into account, fn 
the framework of the stochastic description, that restriction means that the 
system possesses an absorbing and/or reflecting barrier at which the prob- 
ability distribution vanishes. Therefore, the corresponding equations must 
involve boundary conditions. Effects connected with the final size of the sys- 
tem are especially pronounced for small objects and they are encountered 
in many applications of the stochastic processes to systems of high com- 
plexity. Taking into account the absorbing barriers is essential in dealing 
with several population and environmental problems. For example, models 
of species extinction involve the size of a refuge in which life conditions are 
favorable, whereas they are extremely harsh outside [fj. inclusion of size of 
the refuge is also crucial in a model of the infection of liantavirus in the 
deer mouse, based on biological observations in North America [2]. 

Time characteristics of the escaping process involve the first passage 
time density distribution, defined as a probability that the time the particle 
needs to reach the absorbing boundary is within the interval -|- dt), 
provided the particle was initially at a given point xq [3]. The average 
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of that distribution - the mean first passage time (MFPT) - is a useful 
quantity to estimate the speed of transport for systems which are defined 
on the restricted area. The MFPT can be calculated also for the boundless 
systems, if a potential restricts the domain. 

For example, an application of the diffusion problem on the finite inter- 
val to the heat conduction between hot and cold baths is recently of wide 
interest. There are several attempts to link such thermodynamical phe- 
nomena as heat conductivity to the dynamical diffusion. In this context, 
the problem of validity of Fourier's law, as a counterpart to the Fick diffusion 
law, is especially interesting. Applicability of dynamical processes, which 
are characterized by the anomalous diffusion coefficient, became clear when 
the anomalous heat conductivity in classical one-dimensional lattice systems 
has been found [3]. A model, called "dynamical heat channels" I5l[6], can 
be constructed by introduction some simplifications, e.g. by neglecting the 
interactions between particles. Then the dynamics can be handled by the 
decoupled CTRW which implies all kinds of the anomalous diffusion. In 
this particular model, the subdiffusive case requires long tails of the waiting 
time distribution; as a result the average waiting time, as well as MFPT 
[7], is infinite. For the heat conduction process that would mean a perfect 
insulator. 

However, that uncoupled version of CTRW does not take into account 
that, in general, the system may be inhomogeneous, i.e. its parameters 
depend explicitly on the spatial variable. This happens in the complex sys- 
tems where long-range space correlations are important and the medium 
structure is crucial for the system properties. As an example can serve the 
transport on the fractal objects [H [9l |10] and, since fractals are ubiquitous 
in nature, its numerous manifestations in various branches of science. The 
transport coefficients must vary with the position if one describes the dy- 
namical properties of materials containing impurities and defects. Physical 
problems which are considered in this context involve conductivity of amor- 
phous materials, the ionic conductors, dynamics of dislocations, transport 
of a dye in porous materials (quenched disordered media) [TT]. In the case 
of the heat conduction, deviations of model calculations from Fourier's law 
indicate that the asymptotic temperature gradient is nonuniform and they 
point at long-range effects |12j . The MFPT for a process which can cor- 
respond e.g. to the Langevin equation with the multiplicative noise, and 
which is also described by the Fokker-Planck equation with the variable 
coefficient, was calculated in Ref. |13] . 

In the present paper we evaluate the MFPT for a jumping process which 
is a version of the CTRW: it is Markovian and takes into account the spatial 
dependences of the problem by introduction the x-dependent waiting time 
distribution. 
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Effects connected with the finite size of the system are especially pro- 
nounced if the particle performs long jumps, namely for the Levy flights, 
when the second moment of the probability distribution is infinite. If the 
distance between boundaries is small, compared to width of the jump length 
distribution, the tails hardly influence the dynamics and the essence of the 
Levy process remains hidden. On the other hand, presence of the barriers 
makes all the moments convergent. 

The paper is organized as follows. In Sec. II we present the definition and 
main properties of the jumping process. In Sec. Ill the MFPT for the system 
which possesses two absorbing barriers is calculated. The consequences of 
introduction of distributions with long tails (Levy fiights) is analysed in 
Sec. IV and dependence on the process parameters is discussed. The main 
results are summarized in Sec.V. 



2. Description of the process 

The process we consider in this paper is a step-wise one-dimensional 
Markov process defined in terms of the jumping size distribution Q{x) and 
the Poissonian waiting time distribution 

Pp{t) = z.(x)e-^(^)*, (1) 

where I'ix) is the jumping rate [Ijj. The process value x{t) is constant 
between consecutive jumps. Since v depends on x, the process is a general- 
ization of the usual, uncoupled CTRW. The master equation is the following 

= -u{x)p{x, t) + r Q{x - x')u{x')p{x', t)dx'. (2) 

Ot J-oo 

In the following, we assume the scaling form, u[x) = \x\~^ {9 > —1), for 
i'{x) which was applied e.g. to study the diffusion on fractal objects |15| . 
Moreover, it was used to describe the transport of fast electrons in a hot 
plasma [16j and the turbulent two-particle diffusion |17j . 

A natural choice for the distribution Q{x) is the Gaussian: 

Q(x) = ^e--V2.^. (3) 



The corresponding master equation for the jumping process can then be ap- 
proximated - by means of the Kramers- Moyal expansion - by the following 
Fokker-Planck equation [19j 

dp{x,t) ^ d^[\x\-<^p{x,t)] 

dt 2 dx^ ' ^ ^ 
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where a is the width of the distribution Q{x). The solution, with the initial 
condition p{x,0) = 5{x), is given by 



p{x,t) = Cg- 



I |6) /__2j£|f+S ^ 

1^1 ^^Pv a^{2+eyH^ 



(5) 



where Cg 
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The mean squared displacement can be 



directly evaluated: ~ t^+e. Then for 9 S ( — 1,0) the superdiffusion 

emerges, for > we get the subdiffusion. The normal diffusion takes place 
for 9 = 0. Therefore, this Markovian process involves all kinds of diffusion. 

The other form of Q{x) is the Levy distribution which is also stable and 
has the broad, power-law tails (0 < // < 2)[18j. The Kramers- 

Moyal approximation of the master equation ^ produces in this case the 
following fractional equation: 



dp{x, t) 
dt 



df^[\x\ ^p{x,t)] 



(6) 



instead of the Fokker-Planck equation ([4]). The solution of the Eq. ([6]) 
represents the Levy process and it can be expressed in terms of the Fox 
function in the following form [191 [20] 



p{x, t) = -HI'I 



a\x\ 



(l-l/^,l//i),(l/2,l/2) 
(0,1), (1/2, 1/2) 



(7) 



where a ~ The solution ([7|) is correct, and equivalent to the 

solution of the master equation in the diffusion limit of large both x 
and t. Since all moments of order 5 > /i of the distribution p{x,t), Eq.d?]), 
are divergent, the kind of diffusion process cannot be determined from time 
dependence of the second moment. Instead, one can introduce fractional 
moments of the order 5 < fi. Alternatively, the renormalized moment of the 
order /i [19] allows us to characterize the diffusion properties of the system 
in the same way as for the Gaussian case and to distinguish the normal 
diffusion (0 = 0), subdiffusion {9 > 0) and the superdiffusion {9 < 0). 

Presence of the absorbing barriers must modify the probability distri- 
bution ^(x, t) for both choices of Q{x): the distributions dwindle with time 
due to the absorption and the broad tails in the Levy case are cut off. As a 
consequence, all the moments are finite. We discuss those problems in the 
next sections. 
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3. Fokker-Planck equation with the boundary conditions 

We consider a one-dimensional motion wliich is restricted to an interval 
[0, L\. The particle performs jumps defined by the probability distributions 
Pp{t) and Q{x) according to the Eqs.([T]) and ([3|). The end points of the 
interval, and L, are regarded as the absorbing barriers; the probability 
distribution is given by the Fokker-Planck equation ^ with the initial con- 
dition t = 0) = 5{x — xo) (0 < xo < L) and with the following boundary 
conditions 

p(0,t) =p(L,t) =0. (8) 

Eq. ^ for this problem can be solved by separation of the variables. Let 
us assume the particular solution in the form p{x,t) = (j){t)^{x). Inserting 
this ansatz to the Eq.(I3|) yields two equations; the function (p can be easily 
determined: 4'{t) = Cexp(— A^t), where A =const. For the function we 
get the equation 

d^yix) 2A2 . ^ 

+ —x%{x) = 0, (9) 

where y{x) = x~^ip{x). The solution of the Eq.Q can be expressed in terms 
of the Bessel functions Ju{x) in the following form 



y{x) = ^^Jv(.+2) [-^0^/'^'^") m 

which satisfies the condition ^(0) = 0- The second boundary condition, 
i^{L) = 0, allows us to determine the parameter A = A„ by means of the 
zeros of the Bessel function 7„: 

_ + 2) 7n 

2^/2 ^^^> 
The general solution can be obtained by summing up over all values of A^: 

p(X,t}-X Z^^n-'l/(0+2) ^^^(9+2)/2^ j exp { j^{e+2)/2>^ 

(12) 

The form of the constant An follows from the initial condition. The orthog- 
onality property of the Bessel function produces, after some algebra, the 
following expression: 



_ e + 2 VXoJl/(e+2) [ L(e+2)/2 Xq ) 
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The series representation of the distribution p{x,t), Eq. (I12p . is conver- 
gent for all X and t. An example of the time evolution of p{x, t), calculated 
according to the Eq.(|12p in which 40 terms has been taken into account, is 
presented in Fig.l. The distributions shift to the right with time and their 
normalization integral becomes smaller - due to absorption at the boundary 
X = L = 5. 




Fig. 1. The solutions of Eq.Q with boundary conditions ^ for 9 = 1. 

Having the distribution p{x, t) calculated, we can determine the survival 
probability: the probability that the particle is still inside the interval (0, L), 
i.e. it has not yet reach the absorbing barrier. It can be obtained by means 
of the formula S{t) = Jq p{x,t)dx and it determines the first passage time 
density distribution f{t) = —dS{t)/dt. The averaging over that distribution 
produces the MFPT: 

roo pL poo 

T= / tf{t)dt= / dx p{x,t)dt. (14) 
JO Jo Jo 

In the case of our jumping process, the direct evaluation of the integral 
yields 

T f In (^'+2)/2\ 

_ 2x0 ^ '^i/(e+2) [jjB+2j72Xo ) /^_.^1^±^^!^n2 A 

(15) 

where we utilized simple properties of the Bessel function. To obtain the 
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MFPT we need to integrate S over time: 

= i?(9T2r? — [j;,„«(7„)i73 ■ 

Fig. 2 presents the survival probability for some values of 9, both positive 
and negative, calculated from Eq. ljlSp . For the problem without absorbing 
barriers, the case with ^ < corresponds to the superdiffusion, whereas 
with > - to the subdiffusion. The figure shows that the tails are always 
exponential. Moreover, S{t) rises with for any t, as expected. For large 
values of 0, beginning of the curve is flat which means that trajectories can 
hardly escape at short time due to the strong trapping. Since S{t) becomes 
actually exponential, the MFPT is finite for all 6 (see Eq. ()16p ). also for those 
which correspond to the subdiffusion. This result is in contrast to that of the 
decoupled CTRW which predicts divergence of the MFPT in the subdiffusive 
case [7]. More precisely, the decoupled CTRW in the subdiffusive case is 
non-Markovian and it assumes the waiting time distribution in the power 
law form. Then the mean time of a single jump is infinite. For the problem 
with the absorbing barriers, one can derive a formula for MFPT directly 
from the waiting time distribution and the MFPT appears infinite. For the 
process presented in this paper, the waiting time distribution is exponential 
and the subdiffusion results from the x— dependence of its coefficient, i.e. 
from nonhomogeneity of the medium. Introduction of that dependence has 
important physical implications. As regards the application to the heat 
conduction problem, the subdiffusive thermal conductivity becomes possible 
also for the systems which are not perfect thermal insulators and then more 
realistic. 



4. Levy flights between the absorbing barriers 

In this section we analyse the jumping process for the system restricted 
by two absorbing barriers for which the jumping size distribution is given 
by the Levy distribution. We calculate the probability density distributions 
and the MFPT as a function of both the Levy index n and the parameter 

e. 

The MFPT problem for the Levy flights on the bounded domain, both 
with and without a potential, is studied extensively in recent years; be- 
side the MFPT, the first passage time distribution has been evaluated as 
a function of the parameters of the Levy distribution, which, in general, 
can be asymmetric. Since the analytical approach is very difficult in this 
case, most of the studies rely on the Monte Carlo simulations [221 [SSI \2^ . 
Nevertheless, recently an analytical solution to the fractional equation with 
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Fig. 2. The survival probability calculated from Eq. lfTS]) for the values of which 
corresponds to the superdiffusion (0 < 0), normal diffusion {9 — 0), and subdiffu- 
sion (6* > 0). 

the boundary conditions, which describes the Levy flights in a homogeneous 
medium (0 = 0), has been found [25j. 

The Levy distribution represents the general stable distribution and in 
that sense it is a generalization of the Gaussian. It accounts for processes for 
which the second moment of the probability density distribution diverges: 
the standard central limit theorem does not apply in this case. Phenom- 
ena which exhibit distributions with long tails are frequently encountered 
in nature. They are typical for systems of high complexity, in particular 
biological social, and financial ones. Therefore, the theory of the Levy 
flights is widely applicable to problems from various branches of science and 
technology. 

One can expect that presence of the barriers will influence the stochas- 
tic dynamics particularly strong in the case of the Levy processes. If the 
interval length L is small compared to the width parameter a of the jump 
length distribution, the power-law tails of the distribution will not manifest 
themselves. In particular, all moments become finite. 

We assume the jump length distribution in the form 



as well as the Poissonian waiting time distribution ([T]). We determine the 
probability density distribution p[x, t) by means of the Monte Carlo method. 




(17) 
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The Levy-distributed jump-size density has been generated by using the 
algorithm from Ref. [27] • The time evolution of individual trajectories, 
which start with the same initial condition, has been performed by sampling 
consecutive values of the jumping size and the waiting time interval from 
the densities Q{x) and Pp{t), respectively. The final results have been 
obtained by averaging over those individual trajectories. Fig.3 presents the 
time evolution of p{x, t) for both positive and negative 9. Similarly as in the 
case of the Fokker-Planck equation, the distributions terminate abruptly at 
the barrier position and they shrink with time due to the absorption. The 
initial delta function at x = 2 is visible up to a long time. 




Fig. 3. The time evolution of probability density distribution for Q{x) in the form 
of the Levy distribution (fT7|) with the parameters /d = 1.5 and a — 1, calculated 
for two values of 9: 9 = 0.2 (left part) and 9 = —0.1 (right part). The vertical hues 
near the absorbing barrier a; = 4 represent very fast fall of the distributions due to 
absorption. 

In the case of the process presented in Fig.3 the width parameter a of the 
driving distribution Q{x) is large, compared to the interval size L, and the 
results are not sensitive to its tails. On the other hand, in the limit of large 
L = L/a one can expect that p{x,t) converges to the distribution which 
corresponds to the process without absorbing barriers. In this case, the 
resulting distribution should not depend on a and the asymptotics should be 
completely determined by fi. Indeed, Fig.4 demonstrates that for a = 0.003 
the tail approaches the form ~ x~^~'^, before it is cut abruptly at the 
barrier position. For the slightly larger value of this parameter, a = 0.004, 
the power-law asymptotics fails to appear. The tails of p{x, t) become a- 
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dependent for relatively large a because the importance of the tails of Q{x) 
gradually declines with a. 
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Fig. 4. The log-log plot of the tails of probability density distribution p{x, t), cal- 
culated with 9 = 0.2 and /i = 1.5, for two values of the width parameter a. The 
solid line denotes the function ^ x^'^l'^ . 

The same procedure allows us to determine MFPT: the average time T 
the edge of the interval, i.e. either or L, is reached. The dependence of 
MFPT on both /i and a is presented in Fig. 5 for the case = 0. In general, 
T rises with \i because then the probability of large jumps falls. Since that 
effect results from the power-law tails, it is weak (the curves are flat) if only 
the central part of the distribution - similar for all ^ ~ is involved, i.e. for 
large a. On the other hand, small values of g result in a strong dependence 
T{fi), which becomes exponential for a = 0.003. This shape persists for 
even smaller a. 

The parameter is crucial to the speed of the transport; in absence of 
the absorbing barriers, the cases with large values of 9 correspond to the 
slow diffusion. It happens because the traps at large distances becomes 
more effective when 6 increases and the transport is hampered due to the 
long waiting times. Consequently, one can expect that the MFPT will rise 
with 9. This conclusion is illustrated in Fig. 6. Growth of the function T{9) 
is especially pronounced for large values of fi (close to Gaussian case = 2) 
because then the average jump length is relatively small and a large number 
of jumps is needed to reach the barrier. Transition from the negative to 
positive values of 9 - which corresponds to the change of kind of diffusion 
in the problem without absorbing barriers - is smooth for all ^u. The latter 
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Fig. 6. The MFPT as a function of the parameter 9, calculated with a — 0.3, for 
some values of /x. 
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5. Summary and discussion 

We have analysed the one-dimensional, step- wise jumping process which 
is defined on the finite interval, bounded by two absorbing barriers. Since 
the waiting time distribution involves the variable, position-dependent co- 
efficient, the inhomogeneity of the medium has been taken into account. 
That, power-law, dependence is a reason of the anomalous behaviour: the 
diffusion can be both weaker and stronger than normal. The calculated 
MFPT is large for the subdiffusion and rises with the parameter 9, never- 
theless it always assumes a finite value. Physical reason behind the latter 
outcome - which can never be concluded from the traditional, uncoupled 
CTRW models of the subdiffusion is the weakening of trapping with the 
increased distance and, as a result, the effective mean waiting time is finite. 
From the point of view of modelling of the anomalous heat conduction, the 
introduction of the x-dependent waiting time distribution allows us to de- 
scribe the subdiffusive heat transport for realistic systems: for those which 
are not perfect insulators. 

Restrictions imposed on the system by the existence of the absorbing 
barriers are pronounced if we allow for long jumps, i.e. if the jumping size 
distribution Q{x) is of the Levy form. The power law tails of that distribu- 
tion can influence the probability density distribution of the process, p{x, t), 
only if Q{x) is narrow, compared to the distance between the barriers, i.e. 
to the system size. In that limit, the sections of p{x,t) which are close to 
the barrier assume the power law shape in the same form as the tails for 
the problem without the barriers. Conversely, for the broad Q{x) those 
tails are hardly visible and the MFPT is almost independent of the Levy 
parameter /x. The dependence of MFPT on the parameter 9 is similar to 
that for the Gaussian Q{x): large values of 9, for which the transport is 
strongly hampered by the traps, result in large T. Nevertheless, it remains 
finite, for any /x and 9. 

REFERENCES 

[1] C. Escudero, J. Buceta, F. J. de la Rubia, and Katja Lindenberg, 
Phys. Rev. E 69, 021908 (2004). 

[2] G. Abranison and V. M. Kenkre, Phys. Rev. E 66. 011912 (2002). 

[3] H. Riskcn. The Fokker-Planck Equation (Springer- Vcrlag, Berlin, 1996). 

[4] S. Lepri, R. Livi, and A. Politi, Phys. Rev. Lett. 78, 1896 (1997). 

[5] S. Denisov, J. Klafter, and M. Urbakh, Phys. Rev. Lett. 91, 194301 (2003). 

[6] P. Cipriani, S. Denisov, and A. Politi, Phys. Rev. Lett. 94, 244301 (2005). 

[7] S. B. Yuste and Katja Lindenberg, Phys. Rev. E 69, 033101 (2004). 



13 



[8] R. Metzler, W. G. Glockle, and T. F. Nonnenmacher, Physica A 211, 13 
(1994). 

[9] R. Metzler and T. F. Nonnenmacher, J. Phys. A 30, 1089 (1997). 
[10] V. E. Tarasov, Chaos 15, 023102 (2005). 
[11] J.-P. Bouchaud and A. Georges, Phys. Rep. 195, 12 (1990). 
[12] S. Lepri, R. Livi, and A. Politi, Phys. Rep. 377, 1 (2003). 
[13] Kwok Sau Fa and E. K. Lenzi, Phys. Rev. E 71, 012101 (2005). 
[14] A. Kaminska and T. Srokowski, Phys. Rev. E 69, 062103 (2004). 
[15] B. O'Shaughnessy and I. Procaccia, Phys. Rev. Lett. 54, 455 (1985). 
[16] A. A. Vedenov, Rev. Plasma Phys. 3, 229 (1967). 

[17] H. Fujisaka, S. Grossmann, and S. Thomae, Z. Naturforsch. Teil A 40, 867 
(1985). 

[18] V. M. Zolotarev, One-dimensional Stable Distribution (American Mathema- 
tical Society 1986). 

[19] T. Srokowski and A. Kaminska, Phys. Rev. E 74, 021103 (2006). 

[20] T. Srokowski, Phys. Rev. E75, 051105 (2007). 

[21] E. Kamke, Differentialgleichungen, Losungsmethoden und Losungen, Band I 
(Geest&Portig K.-G., Leipzig, 1956). 

[22] B. Dybiec, E. Gudowska-Nowak and P. Hanggi, Phys. Rev. E 73, 046104 

(2006). 

[23] M. Ferraro and L. Zaninetti, Phys. Rev. E 73, 057102 (2006). 

[24] S. L. A. de Queiroz, Phys. Rev. E 71, 016134 (2005). 

[25] A. Zoia, A. Rosso and M. Kardar Phys. Rev. E 76, 021116 (2007). 

[26] B. J. West and W. Dccring, Phys. Rep. 246, 1 (1994). 

[27] J. M. Chambers, C. L. Mallows and B. W. Stuck, J. Am. Stat. Assoc. 82, 704 
(1987) 



